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We introduce a new phase- field model which allows for simulation of incoherent solid/solid trans- 
formations. Contrary to previous models which impose coherency at the interface, the zero shear- 
stress condition characteristic of incoherent solids is recovered in the limit of small interface thick- 
ness. For the sake of clarity, we limit ourselves to the case of stress-driven phase transitions between 
two elastic solids. However, since we use a variational formulation which has a clear thermody- 
namic interpretation, the extension of the model to the cases of chemical and/or thermal diffusion 
is straightforward. 
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In recent years, phase-field models have become in- 
creasingly successful in simulating reahstic microstruc- 
ture formation during phase transitions. These models, 

[ which replace mathematically sharp interfaces by diffuse 
ones provide powerful numerical methods which avoid the 

' tracking of the interface (free-boundary problem). The 
idea can be traced back to studies of phase transitions by 

[ Ginzburg and Landau ^ , Cahn and Hilliard ||] , Halperin 
et al. ^ or Langer and Sekerka ||], but it has received 

[ renewed interest because recent phase-field models make 
quantitative numerical simulations possible. The efforts 

[ of the community in this field can roughly be divided 
into two main directions : studies of solidification and 

' dendritic growth on the one hand, coherent solid/solid 
transformations such as martensitic transformation on 

[ the other hand. Among the recent achievements within 
these two sub-fields one can cite the " thin-interface" limit 

' of Karma and Rappel Q , simulations of well developped 
3D dendrites Pi , studies about grain boundary motion 

[ 0, eutectic |^or peritectic ^ growth and martensitic 
transformation Recently, Karma extended the 

model to include phases with asymmetric transport co- 
efficients jl^ and some authors brought attention onto 
the problem of phase-field modeUing of solid/fluid equi- 
libria with elastic stresses. The latter studies were con- 
cerned with the so-called Asaro-Tiller-Grinfeld instability 
1^ , Q or fracture phenomena . All of this is of course 
only a small sample of the existing litterature aimed at 
demonstrating the great variety of problems which can 

[ be tackled using this method. However, the interested 
reader will be able to find quite a large amount of refer- 

I ences within these papers. 

According to the thermodynamic description of solids 
given by Larche and Cahn [|6|, we can distinguish be- 
tween three kinds of thermodynamic boundary condi- 
tions : solid/fluid, incoherent solid/solid or coherent 
solid/solid. The latter imposes continuity of the trac- 
tions at the interface whereas the other ones impose 
continuity of the normal traction and set the tangential 
component to zero. Only solid/fluid and coherent equi- 



libria have been so far successfully modelled within the 
phase-field approach. Although solid/fiuid and incoher- 
ent solid/solid equilibria appear to have the same kind of 
mechanical boundary conditions, their behavior become 
different when diffusion potentials are taken into account 
Ip^ , [l7| . Moreover, it was noticed by Kassner et al. 
that the direct generalization of their model would yield 
coherency at the interface. Thus, an "incoherent" model 
demands a treatment of its own. 

One of the main applications of such a model would be 
the correct description of intermetallic compounds for- 
mation by bulk isothermal diffusion, a phenomenon for 
which incoherent boundaries are rather the rule than the 
exception. These compounds can sometimes exhibit very 
irregular morphologies, as a consequence of a growth in- 
stability which is hard to explain if elastic stresses are 
not taken into account. This instability has been ob- 
served for example in Ni-Si or Mo-Si binary cou- 
ples. The problem is thus of practical importance since 
silicides are used as conducting layers in MOS (metal- 
oxyde-semiconductor) structures and the necessity to un- 
derstand and control interfacial morphologies is obvious 
in this case. 

In the following, we describe our model and derive its 
sharp-interface limit. This provides us with a Gibbs- 
Thomson condition for equilibrium at the interface. As 
a consequence of this relation, it will appear that a 
" Grinfeld-like" instability is possible at the interface be- 
tween two uniaxially stressed incoherent solids. 

We now introduce the model. Since it is supposed to 
yield quantitative modelling and most thermodynamic 
data are expressed at constant pressure, we choose to 
formulate the model with stresses as primary variables. 
Then, it is easy to derive free-energy contributions un- 
der non-hydrostatic stress conditions, using Maxwell re- 
lations to integrate from a state of constant pressure. 
We consider two solid phases Sa and S/s which obeys 
Hooke's law with elastic coefficients {EonVa) and {Ep^vp) 
respectively. The underlying assumption of linear elas- 
ticity means that the strains have to be small, a con- 



2 




FIG. 1: Evolution of elastic coefficients E and u as a function 
of <j) in the interfacial region G [0,1]. 



dition which we will later make more quantitative. Note 
that this means that the applied stress cr must remain 
much smaller than Young's modulus E. We introduce a 
phase-field (p, along with a stress field aij. We suppose 
that mechanical equilibrium is realized much faster than 
any other thermodynamic process which implies the re- 
lation djaij — 0. In this equation and all subsequent 
developments, Einstein convention is used. Subscripts 
preceded by a comma will denote spatial derivatives. We 
then postulate the following form for the free energy: 



T= dV <T 
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where g{(p) — — 0)^ is the usual double- well potential 
with minima at (j) — (phase Sa) and (p ~ 1 (phase Sp) 
and r — 3j/d, 7 being the surface free-energy, which we 
take to be isotropic. The reference state for strain has 
been chosen as a state of zero stress and strain in both 
phases. The precise form of the functions E{(j)) and ^{(j)) 
will only be needed later in the derivation. We plot them 
schematically in Fig.|l]. Note that we have E{0) — E^, 
E{1) = Ep, v{Q) - v^, = vp, E'{0) = E'{1) = 0, 
i/'(0) = !/'(!) = 0. 

We now have to specify how the fields evolve in time. 
Since we assumed mechanical equilibrium, J-' must be 
stationnary with respect to variations of aij. As cry is 
constrained by the relation djaij = 0, some care is needed 
in writing down the equilibrium condition. Let us define 
the generalized strains Sij: 
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which is just St Venant's compatibility equation express- 
ing the fact that the strain tensor Sij is the symetric part 
of the gradient of a displacement field. 

For the phase-field (f), we impose non-conserved dynam- 
ics (Model A of the classification by Halperin et al. [^). 
We thus have the equation: 



dt 56 
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where R = l/(3A:e) is a phenomenological rate constant. 
Expanding and casting the last equation into nondimen- 
sional form, we obtain: 



dt 



3e 



1 



2E 



2E 



'kk 



Then, it can be shown that the stationnarity condition 



(6) 

with e = djlc and Iq = ^E j is the Griffith's length 
(E average Young's modulus, tr applied stress). Since 
the so-called " sharp-interface limit" , which allows to re- 
cover the usual Gibbs-Thomson equation, is a matched 
asymptotic expansion with e as a small parameter, the 
interface thickness must be much smaller than the Grif- 
fith's length in order for the phase-field model to be use- 
ful. This sets a limit on the range of appHed stresses 
which we can be simulated within this framework. Even 
if other phenomena such as chemical or thermal diffusion 
are present and set a natural length scale, it is always 
possible to find a range of external stresses small enough 
to keep the elastic term in (^ much smaller than the 
double-well potential term, a condition which is neces- 
sary in order to obtain meaningful asymptotics. Indeed, 
the extension to the case of a diffusive transformation is 
straightforward as long as no cross-effects (such as ther- 
mal expansion or compositionally-induced stresses) are 
present: we just need to add the elastic part (||) to free- 
energies used in "diffusive" models. Since there is no bulk 
coupling between stress and temperature or concentra- 
tion fields, they can be treated as independent variables 
and the analysis leading to thin or sharp-interface limits 
is not modified. 

From now on, all length scales will be measured in 
units of /g, which amounts to replacing d by e in (|^). 
Moreover, we will limit the following discussion to the 
two-dimensional case, the three-dimensional case being 
just a straightforward generalization of the latter. 

We now perform the asymptotic expansion. In the 
outer domain, i.e. the domain in which gradients of (p 
vanish, the phase-field equation (^) gives the condition 
(7'(0) = at lowest order in e. This equation yields three 
solutions, = 0, 1 or 1/2, the latter being unstable. 
Since 5', E' and v' vanish when <?!) = or 1, these two 
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solutions are valid at all orders of the expansion. The 
outer equations for the elastic fields then read: 

(1 + Up)(Ti.j^kk + o-kkAj = , p = (a , /3) , (7) 

which are just the usual Beltrami-Mitchell equations for 
each phase. 

Let us now consider the inner equations, which solu- 
tions will allow to match the outer solutions at the in- 
terface. We switch to a curvilinear coordinates system 
comoving with the interface. We choose the r axis to be 
perpendicular to the interface (i.e. parallel to and 
oriented from phase to phase Sa- The other coordi- 
nate s is taken along the tangent to the curve </) = 1/2. 
We then stretch the r variable, introducing p = r/e. Us- 
ing capital letters for the inner fields, we obtain at leading 
order: 

^2^(0) ^ 2g'($(")) (8) 
for the phase-field and: 

apE(°) = , d,i:f} = , a^^i^i? ^ o (9) 

for the elastic fields. From equation (^), we obtain that 
Spp'', Sps^ and i?is'' do not depend on p within the inner 
region, a conclusion identical to the results by Kassner 
et al. 0. 

We now give the detailed form of functions and 
E{(l)). We define the function = 0^(3 - 20). We 

then have : 



E{cl>) ^E^ + {Ep - E^)h{4>) 
1 -I- v{(t,) = {aE{4)) - bf 

where a and b are defined according to: 



(10) 
(11) 



(Ep - E„)a = y^l + iya + ^1 + Vfj 



With this choice of functions, it is easily verified that 
we have v{0) = z/„, v{l) = up, E{0) = E^, E{1) = Ep, 
that the derivatives v' and E' are zero when <j) — Q ov 
1 and that there exists 0c G [0,1] such that 1 -I- v{4)c) — 
v'{<j)c) — 0. Since we have: 



^(0) ^ i + K'i^'"^) 
£;($(o)) 
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and 1 -|- ^{(j)) vanishes somewhere in the inner region, 
r(s, t) and thus Ep? must be zero. We now match the 
inner and outer fields, using the matching condition 



ip{r,s,t)\^-_ 
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valid for any physical field ip. Equation (||) yields the 
classical 'I>'^°-' = (1 — tanhp)/2 solution which obviously 
matches the outer solutions = and 0=1. For the 
stress field, we obtain: 



{Ep -Eo,)b^ Ep^l + j/„ + Ea^^l + vp 



(14) 

which are the mechanical equilibrium conditions for in- 
coherent interfaces. 

At first order in the expansion for $, we obtain an 
equation of the form = y(<I>*^°\ E^^^*) where we 

defined the following self-adjoint operator: L = — 
2g"($(0)). In order for a non-trivial solution to exist, 
the right-hand side must be orthogonal to the zero-mode 
i9^(j)(o)^ This condition is called the "solvability con- 
dition" and reads, using Epg-* = and evaluating 0- 
derivatives at = (JjC^); 



= 



dpdp^^"^ I -39p$("' [kv + 7Kj 



N d 


'l + u' 


) ^ d(j) 


2E 



1 ^^.^ 



2E 



1 -r '^ss 



(15) 



r 



We finally note that is a monotonous function of 
0, which allows for a change of variables in the interval 
G [0,1]. Using equations (||), (^) and ( [Til ) to express 
eIs^ as a function of i?ss(<i>(°)), performing the change of 
variables and computing the integral in (|l5|), we finally 
obtain for the local velocity of the interface v: 




)) +7'« 



(16) 



which is the Gibbs-Thomson equation for our model. In 
doing so, we assumed Epp' =0 for simplicity. This choice 
is natural since we chose the equilibrium state with zero 
stress as the reference state for strains. 

When performing the above analysis in two- 
dimensional space, we implicitly assumed that the elastic 
fields corresponded to a state of plane stress (ctzz = 0). 
The plane strain case [ezz = 0) is easily obtained from 
equation (p6[), replacing E hy E / {\ — v^). In the lat- 
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ter case, we could also assume that phase Sa describes 
a fluid, which gives Ea = 0. Then, we recover exactly 
the same Gibbs-Thomson equation as was obtained by 
Kassner et al. in their study of Grinfeld instability [ p^ . 
This in turn qualitatively shows that a " Grinfeld- like" 
instability is also possible at an incoherent solid/solid in- 
terface. 

In the present state, our model is restricted to pure 
"stress-transitions" whereas an incoherent phase trans- 
formation almost always involves mass transport across 
the interface. This means that the model should at least 
be extended to the case of chemical diffusion, which is 
quite straightforward as explained earlier. Going further 
toward realistic systems, it would be interesting to incor- 
porate stresses induced by compositional inhomogeneities 
into the picture. Such an analysis has been performed 
by Gurtin [^9| for the case of coherent equilibrium, us- 
ing a phase-field model with a conserved order parameter 
(Cahn-Hilliard equation) and a different thermodynamic 
approach. In the present case, the extension is not con- 
ceptually difficult since the variational formulation and 
its thermodynamic interpretation allow for an easy trans- 
lation of Larche and Cahn analysis ||l^ into " phase-field 
language" . However, it can not be decided from such 
qualitative arguments whether a zero-th order expansion 
for the inner elastic fields would be sufficient to recover 
the correct Gibbs-Thomson equation. 

To summarize, we have developped a new phase-field 
model aimed at simulating stressed incoherent solid/solid 
interfaces. In the limit where Young's modulus of one 
of the phase goes to zero, the Gibbs-Thomson equation 
for solid/fluid equilibrium with elastic stresses is recov- 
ered p4) . However, contrary to the model by Kassner 
et al., the present model is formulated with stress as a 
primary field which makes it easier to link to thermody- 
namic data for quantitative simulations. Moreover, such 
a formulation is easier to simulate when applied to sys- 
tems where stress is imposed at external boundaries, the 
Neumann boundary conditions for the Navier equations 
becoming Dirichlet conditions for the Beltrami-Mitchell 
equations. The model lends itself to easy generalization 
to more complex situations (thermal or chemical difu- 



sion, thermo-elastic materials, compositionally-induced 
stresses ...). The specific case of incoherent equilibrium 
with diffusion and chemical self-stresses will be the object 
of a future publication. 
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